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ABSTRACT 

We combine data from the Sloan Digital Sky Survey (SDSS) and the Arecibo Legacy 
Fast ALFA Survey (ALFALFA) to study the cold atomic gas content of galaxies in 
groups and clusters in local universe. A careful cross-matching of galaxies in the SDSS, 
ALFALFA and SDSS group catalogs provides a sample of group galaxies with stellar 
masses IO^'^Mq < M* < and group halo masses < 

1015.o^-1j^q_ Controlling our sample in stellar mass and redshift, we find no significant 
radial variation in the galaxy H I gas-to-stellar mass ratio for the halo mass range in 
our sample. However, the fraction of galaxies detected in ALFALFA declines steadily 
towards the centers of groups with the effect being most prominent in the most massive 
halos. In the outskirts of massive halos a hint of a depressed detection fraction for low 
mass galaxies suggests pre-processing that decreases the H I in these galaxies before 
they fall into massive clusters. We interpret the decline in the ALFALFA detection 
of galaxies in the context of a threshold halo mass for ram pressure stripping for a 
given galaxy stellar mass. The lack of an observable decrease in the galaxy H I gas-to- 
stellar mass ratio with the position of galaxies within groups and clusters highlights the 
difficulty of detecting the impact of environment on the galaxy H I content in a shallow 
H I survey. 

Subject headings: galaxies: clusters: general - galaxies: groups: general - ISM: general 


1. INTRODUCTION 

The amount of cold gas in group and cluster galaxies provides information on the impact 
of these environments on galaxy evolution. Previous studies have shown that the distribution of 
cold gas in galaxies, traced by neutral hydrogen (H I), is affected primarily by tidal interactions 
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(iToomre &: Toomrg|l972l l and ram pressure stripping (jGunn &: Gottlll972l i fsee lBoselli &: Gavazzi 
20061 for review of physical processes in galaxy environment). 


Tidal interactions are expected to be important in galaxy groups but l ess so in clus ters due 


to the large relative velocities of galaxies and short interaction time scales (jMihosI l2004l ) . Alter¬ 


natively, ram pressure stripping is thought to be responsible for H I-deficiency in cluster galaxies 


(|Boselh &L Gavazzil I 2 OO 6 I : iRoedigen l2009l : iGavazzi et al.M2013l b Furthermore, ram pressure strip 


cals 




Freeland fc Wilcots 


Lucero et al. 

2005: 

McGarthv et al. 

20081 and dwarfs (e.g., 

McGonnachie et al. 

2007; 


vation’ (iLarson et al 


201 lli. and the hot halo gas in galaxies fe.g.. IJeltema et al.ll2008l l causing ‘star- 


1980|)- For this work we select our sample to minimize the impact of tidal 


interactions so that we can focus on the effects of ram pressure stripping. 

The impact of ram pressure stripping on galaxies has been well studied in the Virgo and 
Goma clusters. Statistical analyses of the H I gas content of their member galaxies provides ev- 
idence for increasing gas depletion toward th e cluster center for optically-selected galaxies (e.g., 
Havnes et al. 1984 : Havnes k. Giovanelli 1986 : Solanes et al. IH) and for H I-selected late type 


galaxies from ALFALFA ( Gavazzi et al 


201 


The disturbed H I morphologies of individual galax- 


(e.g.. 

Bravo-A 

faro et al. 

2000, 

20Q1 

Ghuno: et al. 

2009, 

2007 

: Growl et ahl 

2005 

Kennev et al. 

2004: 

Vollmer et al. 

20081. Nevertheless. 

Lavlor et al. ( 

2012 

) and 

7ortese et al. 

(20081 find onlv marginal 


evidence for increasing H I-deficiency toward the inner regions of the Virgo and the Goma/Al367 
clusters respectiv ely using H I-selected samples from the Arecibo Galaxy Environment Survey 
(jAuld et al.ll2006l f. This lack of H I depletion in these clusters is likely due to the limited dynamic 
range of this flux-limited H I survey. 

While statistical and individual observations suggest that interaction between galaxies and 
the intracluster medium is important for the remov al of gas from galaxies in the most massive 


groups, the situation is less clear in smaller groups (jRasmussen et al.l 120121 ). In a small number 


of their gas (e.g.. 

Bureau V Garienan 

2003: 

Rasmussen et al. 

2006|: 

Seneupta et al. 

2 OO 7 I 1 


Freeland fc Wilcots 


2011 


McGonnachie et al. 2007 


2007l i. Some of these groups also have X-ray detections of 


a hot intergalactic medium (IGM) consistent with the ram pressure stripping explanation. 

A few recent studies have exami ned the impact of envi ronment for wide range of galaxy group 
halo masses using statistical sample. iGatinella et ahl (l201.3l i used ~ 800 galaxies from the GALEX 


Arecibo SDSS survey fGASS: IGatinella et al 


2010 


201^ ) to show that the H I gas fraction (defined 
as Mhi/M* in their work) is at least 0.4 dex smaller for gal axies with stel l ar ma sses M* > lO^^'^^Mfr 


that reside in groups with halo mass > 10^'^ lEabello et al.1 (|2012l i stacked ALEALEA 


0 


spectra of galaxies meeting the GASS selection criteria in bins of stellar mass and local density and 
found a strong decline in the galaxy H I content in dark matter halos w ith mass > 10^^Mq. 
Hess Sz Wilcots (2013) used the SDSS group catalog ( Berlind et al. 2006) and the ALFALFA 40% 
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catalog (jHavnes et alJ [201l|) to show that as group membership increases, galaxy group centers 
become increasingly deficient of H I-rich galaxies. This work builds on these studies by using a 
control sample (i.e., isolated field galaxies with similar stellar mass and redshift) to compare the 
impact of the groups and clusters to the H I gas content of galaxies. 


Hydrodynamic simulations show that ram pressure stripping can remove the outer H I gas 

from a galaxy falling into cluster (e.g., 

Abadi et al. 

1999; 

Marcolini et al. 

2003 

: Roediger & Briiggen 

to 

o 

o 

Tonnesen & Brvan 

2009, 

20101 within 10-100 Mvr ( 

Abadi et al. 

1999; 

Marcolini et al. 

2003; 


Roedigei 


200S: 


Tecceetal 


2010|) while complete unbinding of the disk gas may ta ke a few 100 Myr 

1 Gvr. lBoselli &: Gavazzi 


dRoedigeilbopg 1. less than the dynamical time scale in groups and clusters ( 


200611 . Theoretical models of ram pressure stripping based on the original formulation in I Gunn &: Gott 


(j 19721 ) find that the ratio of the mass of the infalling galaxy to that of th e host group o r cluster is an 
important parameter governing the efficiency of ram pressure stripping (lHesteiil2006l ) and the frac- 
tion of galaxies wi thout cold gas increases toward the cluster center due to increasing ram pressure 
(jTecce et al.ll2010l ). Therefore one can expect that gas depletion will increase toward the center of 
groups and clusters and that the depletion will be more significant in more massive clusters. 


We use a statistical sample of galaxies with H I measurements from the ALFALFA survey to 
investigate the distribution of galaxy H I gas-to-stellar mass ratio as a function of the projected 
distance from the center of galaxy group halos for wide range of halo masses. We have combined 
the ALFALFA H I source catalog, the SDSS photometric and spectroscopic catalog, and the SDSS 
group catalog. In addition, we have created a control sample to remove the biases that arise 


from the correlati on between stellar mass and H I gas-to-stellar mass ratio (iGatinella et al, 


2012 


Huang et al.ll2012lL 


In Section [21 we introduce the data sets used in this study and discuss our sample selection. In 
Section [3] we describe the procedure used to match SDSS galaxies, group catalog, and H I sources 
from ALFALFA. In Section [4l we examine the distribution of galaxy H I gas-to-stellar mass ratio 
and the fraction of the H I detected galaxies as a function of the projected distance from the group 
center and derive a simple ram pressure stripping criterion depending on group halo mass which 
provides a plausible explanation to our hndings. In Section [5] we discuss the implications of our 
results and the relevant issues and summarize the results in Section [6l If not explicitly stated 
otherwise, we use a spatially flat LGDM cosmology with = 0.3, Da = 0.7 and Hq = lOO/i 
km s“^Mpc“^, where h = 0.73. 


2. DATA 


The galaxies in this sample are selected from the NASA Sloan Atlas catalog El (Blanton et al 


20 111 ) and have H I 21 cm measurements in the ALFALFA 40% catalog (jHavnes et al.ll201ll ) and 


Tttp: //www .nsatlas . org/ 
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a group identification in the SDSS group catalog (|Yang et alJ 120071) . We restrict the sample to 
galaxies in halos with masses more than 4 members, and 

redshifts of 0.01 < z < 0.055. A complementary sample of isolated galaxies is used to create 
a control sample for this study. The cross-matching of the three catalogs is discussed in detail 
in Section [3] but it is useful to note that the catalogs have very different selection functions and 
completeness. Requiring that the galaxies be in all three catalogs produces a sample that is mainly 
limited by the ALFALFA selection function. 


2.1. The Optical Galaxy Sample 


Our optical galaxy sample is selected from the NASA Sloan Atlas catalog (jBlanton et al.ll201ll ) 
which includes photometric and spectrosc opic measurements for ~150,000 galaxies within ~ 220 


Mpc 0.055) selected from SDSS DR8 ( Aihara et al.ll201lh . The size of the NASA Sloan Atlas 


catalog is much smaller than the full SDSS DR8 because of the redshift limit. However, the NASA 
Sloan Atlas catalog provides more reliable and accurate photometry than the SDSS DR8 catalog. In 
particular, the improved backgro und subtraction rem oves the size-dependent flux bias for galaxies 
with half-light radii of rso < lOO^^f Blanton et al.ll201lh . The NASA Sloan Atlas catalog also provides 
an improved photometric center for galaxies in which the SDSS spectrum is not coincident with 
the photometric center and has eliminated objects where the SDSS spectroscopic object is a sub¬ 
clump of a larger galaxy. These improvements allow for a more reliable estimate of stellar mass 
and improved matching with the ALFALFA catalog galaxies. 

The optical photometry from the NASA Sloan Atlas catalog is used to derive the stellar 
masses of t he ga laxies. Galaxy stellar mass is computed from galaxy color and luminosity following 
Bell et ahl 


log 


M* 


= -0.306 1.097(5 -r)- 0.1 - 0.4(M,. - 51og/i - 4.64) 


( 1 ) 


h-^MQ 

where Mj. is the SDSS r-band absolute magnitude and {g — r) is the g — r color based on the Sersic 
model magnitude. The color-based stellar mass measurement is used for consistency with the stellar 
mass definition in the SDSS group catalog (jYang et al.l 120071 ). However, the results presented here 
do not change if stellar masses derived from SED fitting are used instead. 


2.2. The H I Galaxy Sample 


The ALFALFA survey is an extragalactic H I 21 cm line su rvey that provides info rmation about 
the neutral atomic gas content of nearby galaxies {z < 0.06, iGiovanelli et al.ll2005l j. The survey 
covers ~7000 square degrees of sky wit h 10 km s~^ velocity re solution (after Hanning smoothing), 
an rms sensitivity of 1.8 mJy beam“^ ( Giovanelli et al.ll2005l ). and a positional accuracy of ~24'' 
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for sources with S/N > 6.5 (|Giovanelli et al.ll2007l h The curren tly available source catalog, a.40, 
contains ~ 15000 sources within 40% of the total survey area (jHavnes et al.ll201lh and includes 
measurements of redshift, H I fl ux, and H I line wid th that can be used to determine the H I mass 
of the source using the relation (jHavnes et al.ll201ll h 


M„t = 2.36 X / Sdv Mp 


0 


( 2 ) 


where D is the distance to the galaxy in Mpc, S is the H I 21 cm flux density, and f Sdv is the H 1 
flux of the source integrated over the velocity. 

Both “code 1” sources, which are highly reliable detections with S/N > 6.5, and “code 2” 
sources, which are lower signal-to-noise detections (4.5 < S/N < 6.5) that are likely to be rea l 
because there is a known optical counterpart at the redshift of the H I source (jHavnes et al.ll201lh . 
are included in our sample. The majority (76%) of the sample consists of “code 1” sources and 
the fraction of “code 2” sources does not vary significantly with the projected distance from group 
center. 

Figured] shows the relationship between the velocity integrated H 1 flux and the H 1 21 cm line 
width measured at 50% of the peak flux , for the sample in this study. The 25%, 50%, and 90% 
completeness limits for the a.40 catalog ( Havnes et al.ll201ll ) are indicated by the red lines in the 
figure. The blue triangles are the group galaxies and the grey points are the isolated galaxies used 
as controls (see TS.ll for details on how this sample was selected from the cross-match of the three 
catalogs). The fraction of the sample below each completeness limit is nearly same for the isolated 
and group samples: 41% and 6% of them are below the 90% and 25% completeness limit. 


2.3. The Galaxy Group Catalog 


The goal of this study is to understand the impact of the group environment on the H I gas 
content in galaxies. To determine the group membership of gal axies in our samp le, we use the 
SDSS DR7 group catalog updated from the DR4 group catalog ( Yang et al. 2007). The catalog 
uses galaxies in the SDSS DR7 spectroscopic sample with 0.01 < z < 0.2 and redshift completeness 
C > 0.7. The group finding algorithm has been extensively tested using mock galaxy redshift 
surve y catalogs and ha s proven to be successful in associating galaxies that reside in a common 
halo (jYang et al.M2005l ). In particular, this halo-based group finder works well for poor groups 
and identifies groups with only one member (i.e., isolated galaxies). The group halo masses are 
determined down to M/j = using two methods: ranking by luminosity and stellar 

mass of galaxies in groups. Although we used the luminosity ranked group halo mass, the results 
do not change if the stellar mass ranked halo mass is used instead. The group finder has been 


^http://gax.shao.ac.cn/data/Group.html 
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show n to correctly select more than 90% of the true halos with > W^^h~^MQ (|Yang et al 
20071 1. which allows us to reliably study our galaxy samples within groups and clusters with halo 
mass lO^^'®/i“^M 0 < < W^^h~^MQ. 


In order to investigate galaxy H I properties as a function of the projected distance from the 
group center, we normalize each galaxy’s radial position by the virial radius of the group. For virial 
radius of group with halo mass M^, we adopt t he radius Riso that encloses an overdensity of 180 
relative to the critical density (|Yang et al.l 120071 1: 


Riso = 1.26/i ^Mpc 


Mh 


10i4/i-iMo 


1/3 


(1 + ^group) 


-1 


(3) 


which is based on the WMAP3 cosmological model parameters (ISpergel et al.l 120071 ). = 0.238, 


Ha = 0.762 and Hq = lOO/i km s“^Mpc~^, where h = 0.73. While these parameters differ slightly 
from those used in this study, using this formulation does not make a significant difference given 
the low redshifts of our sample (0.01 < z < 0.055). 

Figure [2] shows the normalized positions of the galaxies in the most massive and least massive 
group halos (M^ < h~^M q and Mh > h~^M q) in this study. Complete details re¬ 

garding the galaxy samples in groups are discussed in Section 13.21 Compared with the distribution 
of the optically-selected galaxies (grey dots), the ALFALFA detected galaxies (blue triangles) are 
rarer near the centers of the groups with the effect significantly more pronounced fo r the more mas¬ 
sive g roups. This lack of H I detections near group centers matches the findings of lHess &: Wilcots 
( 2013l L In Section [4.21 we investigate the H I detection fraction more carefully by using a control 
sample to avoid the effect of H I selection bias by the masses and distances of galaxies at different 
positions within group halos. 


3. CROSS-MATCHING AND SAMPLE DEFINITION 
3.1. Cross-matching SDSS, ALFALFA, and the Group Catalog 

For each SDSS galaxy in the NASA Sloan Atlas catalog, the associated group halo was identified 
by locating the SDSS galaxy in the group catalog that is within 5" in position. This produces 
103,287 SDSS galaxy matches. Since galaxies in the SDSS group catalog and in the NASA Sloan 
Atlas catalog are the same sources, 5" positional matching is very robust and the redshift difference 
is typically Az ~ 0.0001 (i.e., 30 kms“^ in velocity). Only 0.1% of matches have redshift differences 
of 0.0003 < Az < 0.0009 (none have Az > 0.0009). The larger velocity differences occur for galaxies 
where the NASA Sloan Atlas Catalog uses a non-SDSS redshift measurement that is deemed more 
reliable. Generally this happens when SDSS does not have a spectrum near the photometric center 


^All matching was done using the catalog manipulation software Topcat (ITavlodLOOSh 



















-7 - 


of the galaxy. All 103,287 NASA Sloan Atlas galaxies that are associated with groups were then 
cross-matched with the a.40 catalog. Sources within 40" in position and 0.001 in redshift (300 
kms“^ in velocity) were considered potential matches. 


In total there were 6444 potential matches between SDSS group galaxies and ALFALFA 
sources. However, since the Arecibo beam is large (~3.5'diameter), there can be more than one 
galaxy in the beam that will influence the measurement of the H I mass. Figure [3] shows an extreme 
example in which there are 4 galaxies in Hickson Compact Group 69 within the ALFALFA beam. 
In this case, the H I measurement is likely not that of a single galaxy. Cases of galaxy confusion 
will be more common in groups than in the field and thus the H I mass of the best matched galaxy 
among the multiple galaxies within the Arecibo beam will be biased (e.g., the best match indicated 
by the cross in Figure [3]). To address this problem we exclude from the matched sample 520 galax¬ 
ies that have additional SDSS sources within a 1.75' radius. Removing these galaxies with close 
neighbors decreases our sample size to 5924 (a decrease of 8%) and removes systems that are likely 
to be experiencing strong and complex gravitat ional interactions as is the case in compact groups 
(e.g.. ICluver et al.ll2ni3l : iRasmussen et al.ll2ni2l ) rather than purely galaxy-IGM interactions. This 
proximity limit excludes galaxies that are within 20 kpc at = 0.01 and 105 kpc at z = 0.055. 


The sample selected from matching the three different catalogs, excluding galaxies with close 
neighbors includes 5924 galaxies. For these galaxies, the catalogs provide SDSS galaxy parameters, 
SDSS group parameters, and the ALFALFA H I 21 cm parameters. The lower limit on the redshift 
{z = 0.01) of the cross-matched sample is set by the SDSS group catalog and the upper limit 
{z = 0.055) is set by the NASA Sloan Atlas catalog. This is the matched sample catalog from 
which we construct the final group sample with controls as discussed in the following section. 


3.2. Defining the Group and Control Samples 


Our cross-matched galaxy sample is limited primarily by the ALFALFA detection limits thus 
biasing the sample towards gas-rich systems. In particular, the detectability of an H I source 
depends on the g alaxy redshift arid the H I r nass which is , at le ast loosely, correlated with the 
stellar mass (e.g., ICatinella et al.l l201fll . l2012l : iHuang et al.l l2012l ). To mitigate these biases, we 
constructed a control sample of isolated field galaxies. The control galaxies were selected to have 
similar stellar masses (AlogM* < 0.1) and redshifts {Az < 0.001) as their group counterparts. 

The catalog of 5924 galaxies provides the sample from which we drew both our isolated and 
group galaxies. Group galaxies were selected to have halo masses larger than 10^^'^/i_^^ALQj_sig¬ 


nifica ntly larger than the 90% completeness limit (lO^^'^/i ^Mq) in the group catalog (jYang et al 


20071 ). and group membership of more than 4 galaxies. Isolated galaxies were identified as galaxies 
residing in groups with one member. The number of galaxies with one group member (i.e., isolated 
field galaxies) and with group membership of more than 4 galaxies are 4372 and 618 respectively. 
The remaining 934 galaxies which are not included in this study reside in groups with 2-4 members. 




















H I gas-to-stellar mass ratio is known to be correlated with g alaxy stellar mass down to 
= IO^^Mq with a large scatter (0.43 dex, ICatinella et al.N2012l b This trend may continue 


to the lower masses but is harder to quantify because of sample incompleteness (jCortese et al, 


201 ll : iHuang et ahl l2012l : Rosenberg et al. in preparation). Because of the large scatter in this 


relationship, a small control sample will introduce a sampling bias and varying the numbers of 
controls will produce differences in the sample variance. Therefore, we use the same number of 
control galaxies, 8, for every group galaxy. 


Forcing each group galaxy to have the same number of controls means that for a larger number 
of controls, the number of group galaxies with enough controls that meet the selection criteria 
(AlogM* < 0.1 and Az < 0.001) decreases. Alternatively, decreasing the number of control 
galaxies increases the number of group galaxies with a full complement of controls, but it also 
increases the variance in the control sample. To optimize the control sample size, we investigated 
how the standard deviation of the H I gas-to-stellar mass ratio changes as the size of the control 
sample changes. For 50 randomly selected pairs of values for stellar mass and redshift that lie within 
the survey range a control sample was assembled. The average log(MHi/M*) for each group of 
controls was calculated, resulting in 50 values of this parameter. The mean and standard deviation 
of these 50 values was then computed. This procedure was repeated for control samples with 
N = 3,5, 7, 15, 25, and 30 and the values for the standard deviation were plotted as a function of 
the number of controls in Figured! The same computation was carried out 20 times and the results 
are shown by the 20 black lines. The red line is the average of the 20 trials. This figure shows that 
small control samples deviate significantly from the expected statistical behavior of variance with 
sample size (shown by the thin solid line), n, which is proportional to Xj^pn, because the controls 
are not representative of the H I gas-to-stellar mass ratio distribution. Alternatively, as the number 
of controls increases, the number of available group galaxies declines. 


Figured! indicates that the optimal control sample size where the variance goes as XjpH while 
maximizing the number of group galaxies is between 7 and 10. We use 8 control galaxies for 
each group galaxy as a balance between maximizing the number of controls and maximizing the 
number of group galaxies. Group galaxies with fewer than 8 controls meeting the selection criteria 
(AM* < 0.1 dex and Az < 0.001) are not included in the final sample. The best controls were 
selected by ranking the isolated galaxies based on the normalized distance from the group galaxy 
in the stellar mass-redshift plane (i.e., dnorm = v^(AlogM*/0.1)2 -|- (Az/O.OOl)^). The 8 closest 
isolated galaxie s were selected as the control sample for each group galaxy. This technique is similar 
to that used by Ellison et al. (2008). 


Figure O shows the distribution of stellar masses and redshifts for galaxies in the group and 
isolated samples. Small dots are the isolated galaxies used as controls and blue circles are the 
group galaxies. The filled circles show the final galaxy sample (390 group galaxies that each have 
8 controls). Note that there are not many group galaxies with a full complement of controls above 
Af* = 10^*^'®^ Mq. For this reason, our hnal sample is limited by M* < 10^®'®^ Mq. The stellar mass 
and redshift distributions of the galaxy group sample and the sample of their associated controls 
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are drawn from the same population wi th more than 99% and 95% confidence respectively, based 
on the Kolmogorov-Smirnov (K-S) test ([Press et alJll986l i. 


Figure [6] and [7] show galaxy stellar mass and redshift for the group and control samples as 
a function of the projected distance from the group center normalized by the group virial radius 
(Riso)- The three panels show different halo mass bins. The AlogM* values are the difference 
between the log of the stellar mass of the group galaxy and the log of the average stellar mass 
of its 8 controls. The Az values are the difference between the redshift of the group galaxy and 
the average redshift of its 8 controls. Error bars associated with each group galaxy correspond to 
the range of stellar mass and redshift spanned by the galaxy’s 8 control samples. The black solid 
lines connect the average of AlogM* and Az in 6 radial bins up to 1.2 Riso with errors defined by 
standard deviations of AlogM* and Az in each radial bin. The dashed lines show the 5*^ and 95*^ 
quantile of the range of radial variation of the average values of AlogM* and Az. This range is 
obtained by randomly shuffling the galaxies’ projected distance within the full range of projected 
distances of the group sample while keeping the values of A log M* and Az the same. These ranges 
provide a measure of the statistical significance of the observed radial variation. Note that the 
error bars for each galaxy are always smaller than 0.1 in AlogM* and 0.001 in Az since the control 
galaxies are selected to have smaller stellar mass and redshift difference than these values. 


The radial variation of the average stellar mass compared to the control sample is within 0.01 
dex which corresponds to less than 3% difference in stellar mass. The radial variation of redshift 
difference is within 0.0002. The stellar masses and redshifts of the group galaxies do not show any 
systematic radial variation with respect to the control sample, which ensures that any measured 
systematic trends as a function of projected group centric distance are not driven by the radial 
variation of these properties. 


4. THE INFLUENCE OF THE GROUP ENVIRONMENT ON THE H I 

CONTENT IN GALAXIES 

The goal of this work is to study the impact of environment on the H I gas content of galaxies 
in groups and clusters. To examine the impact of environment, we compare the gas content of 
galaxies within groups and clusters to that in field galaxies (i.e., the control sample) as a function 
of projected distance from the group center. 


4.1. Radial Distribution of H I Gas-to-Stellar Mass Ratio as a Function of Halo 

Mass 

Figure [8] shows Alog(MHi/M*) as a function of projected distance from the group center nor¬ 
malized by i?i8o- The value Alog(MHi/M*) is defined as: 
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Alog(MH:/M*) = log(MHi/M*) 

group log < (Mhi/M*) 

control ^ ■ (4) 

where < control > is the average value of H I gas-to-stellar mass ratio of the control 

sample. 

The three panels show the results for three different halo mass bins, Mq< < 

1013.35/^-1 1013-35/1-1 Mq< Mh < 1013-85/j-i Mq and 1013-85/1-1 Mq< Mh < IO 13 -O/ 1-1 Mq. 

The halo mass bins are chosen so that there are similar number of galaxies (~ 120 — 150) in each bin 
and so that the variation in the number of galaxies across the radial bins is small (i.e., the standard 
deviation of the number of galaxies in the 6 radial bins is ~ 50 — 58% of the average number of 
galaxies). The results presented here do not change significantly as the boundaries between the 
bins are shifted from 10i3-3*ih“i Mq to 10i3-'i*ih“i Mq and from 10i3'83/i-i Mq to 10i3'®i^/i“i Mq. 
The error bars and symbols in this hgure are the same as in Figure [6l The black points connected 
by a solid line show the average Alog(MHi/M*) in each bin of projected distance from the group 
center. 

Overall there is no statistically signihcant variation in the average A log(MHi/M*) as a function 
of the projected group centric distance. However, there is a hint of a decreasing Alog(MHi/M*) 
for galaxies in the central region (d/Riso < 0.2) of the group halos in the most massive halo bin 
0. The average Alog(MHi/M*) in the smallest projected distance bin is slightly smaller than the 
range enclosed by the 5*^ and 95*^ quantile of the variation but the deviation is small so it can only 
hint at the impact of the most massive halos on the gas content of galaxies near their centers. 


4.2. Radial Distribution of Detection Fraction as a Function of Halo Mass 

The lack of a statistically significant result in Figure [8] may be driven by the selection of 
only the most gas-rich galaxies by ALFALFA. Therefore, we investigate the fraction of galaxies 
detected by ALFALFA as function of projected distance from the group center. For this analysis 
we compute the detection fraction for the group galaxies and compare it with the detection fraction 
of the associated controls. In practice, we define “detection ratio” as the ratio of the fraction of 
group galaxies detected at a given distance from the group center to the fraction of the associated 
control galaxies that are detected: 


Detection Ratio = ^^det/Ntot)group 

( ^det / ^tot ) control 


(5) 


^These massive halos were identified with some of known galaxy clusters in the Abell catalog ( Abell et al. 
including A160, A996, A1185, A1656 (in Coma), A1795, A1797, A2040, A2052, A2063, A2147 (in Hercules) and 


1989li . 


A2592. 
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where A^^et is the number of galaxies detected by the ALFALFA survey and Ntot is the total number 
of SDSS galaxies in the ALFALFA survey region. 

For each bin of projected group centric distance, we compute the numerator {Ndet/^tot)group 
using all group galaxies that meet the selection criteria dehned in Section [3.21 Then, for each SDSS 
group galaxy (both detected and not detected by ALFALFA) used for computing the numerator 
{^det/^tot)group-, we identify a set of 8 isolated control galaxies (number of group members equals 
one in the group catalog) with similar stellar mass and redshift (AM* < 0.1 dex and Az < 0.001). 
The denominator in Equation [5] is the fraction of the ALFALFA detected control galaxies in the 
projected group centric distance bin. The denominator provides a normalization of the detection 
fraction that is sensitive to the mass and distance of the group galaxies in each bin. For the control 
sample, d/i^igo and halo mass do not have a physical meaning. However, each control galaxy is 
associated with an individual galaxy that does have a projected separation from the group center 
so each control is associated with a particular projected separation through its associated galaxy. 

different ways to select a control sample (e.g., the single best match like 
or allowing control galaxies to be selected only once), we performed this 
analysis using different selection criteria for the control sample and find no significant difference 
in the results. Therefore, we present the results for 8 controls, consistent with the analysis from 
Section 14.11 


Because there are 


m 


Ellison et al.l (201a) 


The upper left-hand panel of EigureO shows the detection ratio for galaxies in group halos as 
a function of the group-centric distance for the same three halo mass bins as in Eigure [HI Error 
bars are Poisson errors. The average detection fraction for the control galaxies is 27% and there 
is little variation among the three halo mass sub-samples (as expected since these are all isolated 
galaxies and are only associated with a halo mass bin because of the galaxies for which they serve 
as a control). When broken down further, there is a 7% variation in the detection fraction of the 
control samples across the radial bins with which they are associated (the higher variance is due 
to the smaller number of control galaxies associated with each radial bin). 

Eor all halo masses there is a decrease in the detection ratio closer to the center of the groups. 
The decrease in detection ratio is signihcantly stronger for galaxies in the most massive halo bin 
reaching only 6% of the detection fraction of the control sample in the smallest projected distance 
bin at d/Riso = 0.1. 


We can only measure the projected distance of galaxies with respect to the group center and 
so in each projected radius bin there are foreground or background “interlopers”, group galaxies 
at significantly larger real separations from the group center than their projected separations. 
Therefore, the detection fraction in a given projected separation bin is an upper limit because these 
interlopers are further from the group center and thus more similar to the field galaxy population. 
To estimate the popula tion of these interlopers in each bin , we use the projecte d mass profile from 


Lokas fc MamonI ([20011) which is based on an NFW profile ([Navarro et al.lll997[ ) with concentration 


parameter c = 7. At projected distances of d/Riso = 0.1,0.3, 0.5,0.7, 0.9, and 1.1, these model 
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predicts interloper fractions of 0.33, 0.16, 0.11, 0.07, 0.06, and 0.05 respectively. 

In the inner-most radial bin d/iiigo = 0.1, the numbers of group galaxies detected by ALFALFA 
are 16, 25, 9 and the total number of group galaxies are 362, 366, 486, for the lowest, middle, and 
highest halo mass bins respectively, which corresponds to detection fractions of 4%, 7% and 2% 
respectively. If one third of galaxies in the inner-most radial bin are interlopers (120, 120 and 161 
for the each halo mass bin), then even a modestly higher detection rate for these galaxies that only 
appear to be in the group because of projection effects would account for many of the detections. 
More relevant than the specific numbers, which are highly affected by small number statistics, is 
that these low detection rates in the centers of groups and clusters are likely to be upper limits 
because of these interlopers. In addition, the sample has been selected to avoid the densest groups 
(no more than one galaxy in the beam of the radio telescope) where some of these effects are likely 
to be even more severe. 


4-2.1. Effect of Galaxy Stellar Mass 

The differences in the depth of the gravitational potential well mean that the impact of envi¬ 
ronment on the H I in a galaxy may be correlated with galaxy mass. The three additional panels 
of Figure [9] show the detection ratio broken down by the galaxy stellar mass for each halo mass 
bin. The galaxies are divided into high (larger red symbol) and low (smaller blue symbol) stellar 
mass systems at a mass of 10 ®'®Mq which roughly corresponds to the median of the galaxy stellar 
mass distribution in the group sample. The detection ratio for low stellar mass galaxies (~ 0.25) 
is significantly less than the ratio for high stellar mass galaxies (~ 0.55) inside the virial radius 
{d/RisQ < 0.8) for the lowest halo mass bin (upper right-hand panel of Figure [9]). The impact of 
the stellar mass of the galaxy on the detection ratio becomes less significant with increasing halo 
mass. For the highest mass halos, the detection ratio for low mass galaxies (M* < 10^'^Mq) just 
beyond the virial radius is only 43% of that in the field, corresponding to a 2a deviation. Although 
the statistical significance is not large, this result may indicate that small galaxies in the outskirts 
of massive halos are already gas poor compared to the field galaxies with similar stellar mass and 
redshift. 


4-2.2. Effect of Galaxy Type: Sersic index 


The morphology-density relation (lDressleiill98rill implies th at there are more early-type galaxies 
towards the centers of galaxy clusters (e.g., Blanton et al.l[2nn5 ). To examine whether the decreasing 
detection ratio in Figure [9] is due to decreasing amounts of gas in the galaxies or to different galaxy 
morphologies as a function of position in the group the three additional panels of Figure [lO] show 
the detection ratio broken down by the galaxy Sersic index for each halo mass bin. The galaxies 
are divided into high (larger red symbol) and low (smaller blue symbol) Sersic index systems at 
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n = 2.0 and then control sample for each galaxy with n > 2.0 (and n < 2.0) is drawn from the 
pool of isolated galaxies with n > 2.0 (and n < 2.0), using the same stellar mass and redshift 
criteria that we us ed. This cut in Sersic index provides a very clean sample of late-type galaxies 
([Mailer et al.ll2009l ) that can be examined for gas depletion. However, the n > 2 sample is more 
mixed and thus more difficult to interpret. We also controlled Sersic index when making control 
sample by restricting Sersic index difference between group galaxy and controls (i.e.. An = 2), 
however found that the result has not been changed qualitatively. The detection ratio for both low 
and high Sersic index galaxies decreases toward the group center indicating that even in the disk 
galaxies gas is being depleted towards the center of the groups. 


4.3. Halo Mass Dependent Ram Pressure Stripping 


The strongly decreasing detection ratio of galaxies towards the center of the most massive 
halos compared with that of the lower mass halos suggests a gas depletion process that is depen¬ 
dent on halo mass. The H I gas depletion seen in galaxies that re s ide in massive clust e rs is best 
explained by ram pressu re stripping (e.g., iBoselli fc Gavazzil 120061 : ICortese et al.N2nill : iRoediger 
20091 : IVollmer et al.l l200llb Some galaxies in gr o ups also represent strong c andidates for ram pres¬ 


sure stripping ('e.g.. iBureau fc CarignanI 120021 : iMcConnachie et al.l 120071 : iR.asmnssen et al.l 12006 


Sengnpta et al.l 120071 ) although it is harder to distinguish t he impact of ram pressure stripping 


from that of tidal interactions in these lower mass halos (^e.g.. iKilborn et al.ll2009l : iRasmussen et al 


20081 1. The decrease in the detection fraction for galaxies with small {n < 2) SctsIc indices implies 
that gas is being removed from the system regardless of the galaxy type and that mergers that can 
transform galaxies from late to early types are not the dominant reason for this result. This gas 
removal can take two forms in these environments, ram pressure stripping (removal of cold gas from 
the disk of the galaxy) and starvation (removal of hot gas from the galaxy halo). Since starvation 
can not be probed directly with these H I observations, we focus our discus sion on ram p r essure 
stripping. The simple criterion for ram pressure stripping hrst proposed by iGnnn fc Gotti (|l972l ) 


is: 


PicM^^ ^ 27rG5]*SigM (6) 

where picM is the density of the hot intracluster medium (ICM), v is the speed of infalling galaxy, 
S* is the galaxy stellar surface density, and Sism is the galaxy interstellar medium (ISM) surface 
density. This criterion implies that gas stripping occurs if the ram pressure (i.e., force per unit 
area) is greater than the gravitational binding force per unit area in the ISM. 

We provide a simple analytic calculation based on the Gunn and Gott (GG) criterion to 
examine the effect of halo mass on ram pressure stripping. The GG criteria has previously been 
modified to be expressed in terms of specific observables: for example, IGM electron density, a 
component of the galaxy velocity relative to the group, central stellar velocity dispersion of the 
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galaxy, and H I gas density of the galaxy (e.g., iFreeland &: Wilcotsll201ll : iGrcevich &: Putmanll2009l : 
Grebel et alJl2003l l. To study the impact of halo mass on ram pressure stripping for a range of 
galaxy stellar masses we use the virial theorem and galaxy scaling relations to rewrite the left-hand 
side of the GG criterion to have explicit dependence on group halo mass and IGM density and the 
right-hand side to depend only on galaxy stellar mass (see Appendix A for details). With these 
changes the GG criterion depends on the mean electron number density in groups or clusters, the 
halo mass and the galaxy stellar mass: 


n 


10 ^cm ^ 


Mh 


1014/i-iMo 


2/3 


> 1.24 


M* 


(1OIOM0 


0.58 


1 -F0.25 X 


M* 


IQiOMo 


(7) 


Massive halos have high IGM densities (e.g., iGavagnolo et al.ll2009l ) due to their strong grav¬ 
itational potentials while small mass halos tend to have low IGM densities. The hal o mass is also 
relat ed with the X-ray luminosity through an observational scaling relation (e.g., lArnaud et al. 
20071 1. Therefore, the virial theorem and X-ray scaling relation imply that the mean electron num¬ 
ber density of the IGM, the halo mass, and the kinetic temperature of the IGM are related by (see 
Appendix A for details): 


n 


10 ^cm ^ 


0.68h'^ 


keV 


-1 


Mh 


1014/i-iMo 


0.824 


( 8 ) 


By combining these two relations and neglecting the redshift dependence, one can determine 
where the GG criterion is valid in the group halo mass and galaxy stellar mass plane for a given 
IGM temperature. Figure [11] shows the distribution of the group sample (black points) in the 
halo mass and galaxy stellar mass plane for galaxies within the virial radius of their group halo 
(where ram pressure is expected to play a role) with black lines representing Equation [3 for IGM 
temperatures ranging from 1 keV to 8 keV. The colored regions show the average Alog(MHi/M*) 
of galaxies in each grid region weighted by their errors in Alog(MHi/M*). Regions with 0-1 galaxy 
have not been colored. The region above the black line for a given IGM temperature is where ram 
pressure wins over the surface gravity of the galaxy. Galaxies in halos with Mh > 
are mostly located in the region where ram pressure wins, for a range of typical IGM temperatures 
(1-8 keV). Alternatively, galaxies in lower mass halos generally reside in the region where galaxy 
surface gravity wins over ram pressure. The galaxies with the lowest Alog(MHi/M*) reside in the 
two bins where ram pressure stripping is expected to play the largest role. 


Figure [TT] implies that for the stellar mass range of the galaxies in our sample W^'^Mq < 
< IO^^'^Mq, and typical IGM temperature (1-8 keV), ram pressure is efficient for halo masses 
larger than ~ which explains the decreasing detection ratio toward the group cen¬ 

ter pr ominently seen in the most massive halo bin Mh > 10 ^^'®^/i“^Mq. Similarly iGatinella et al 


( 2013 1 find that massive galaxies (M* > 10^*^ Mq) in halos with Mh > 10^^ Mq show a decre¬ 
ment of ^0.4 dex in H I gas-to-stellar mass ratio relative to the similar galaxies in lower density 
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environments. Using semi-analytic galaxy formation models, iFabello et al.l (|2012l i also infer that 
H I-deficient galaxies are more common in dark matter halos with > 10^^ The modified ram 
pressure stripping criterion that depends on halo mass is used here to show that ram pressure strip¬ 
ping can be effective in depleting the H I for wider range of stellar mass IO^'^Mq < M* < 


DISCUSSION 


5.1. Gas Removal in Group and Cluster Halos 

Using the SDSS and ALFALFA surveys, we have examined the evidence for ram pressure 
stripping of H I gas in groups and clusters. We find no statistically significant evidence for a 
change in the galaxy H I gas-to-stellar mass ratio as a function of projected dist ance from the 


group center in this sample desp i te the existing evide nce for gas depletion in clusters (jHavnes et al 


19841 : iHavnes &: Giovanellil 119861 : ISolanes et al.l [200l|) . This result i s, however, consist ent with the 
previous H I studies of gas depl etion in Virgo Clus ter galaxies (|Tavlor et al.l l2012l i and Coma 
Supercluster galaxies (Abell 1367. ICortese et al.ll2008l i. The distribution of H 1 gas-to-stellar mass 
ratio in galaxies shows a large scatter (e.g., 0.43 dex in ICatinella et al.l 120121 ) even at fixed stellar 
mass which drives a large standard deviation (approximately 0.4 dex or a factor of 2.5 as seen 
in Figure 11]) in Alog(MHi/M*). Although, in principle, it is possible to detect a difference in the 
average H I gas-to-stellar mass ratio between groups and field with enough statistics, our sample 
of groups with 8 controls (390) is not enough to detect the difference. A robust detection of gas 
depletion in groups will require deeper H I surveys that increase the dynamic range of detected gas 
fractions for a large number of galaxies in groups and clusters. 

While the direct measurement of gas removal in groups and clusters using ALFALFA detected 
galaxies alone is not possible without relyin g on deeper H I survey for a sam ple with stellar mas s 


threshold fGASS: ICatinella et al. 


2010 , 


20]j ) or stacking of ALFALFA spectra ( Fabello et al.ll2012l ) 


we confirm the results of lHess fc Wilcota (j2013l ) using a sample with matched controls and measure 
this removal indirectly through the decreasing detection ratio, which is a ratio of detection fraction 
for group galaxies to that for isolated galaxies. By controlling our sample for stellar mass and 
redshift we show that this is a real effect not a function of selection effects and sampling. We find 
that most significant decrease in detection ratio occurs in the most massive halos, as expected from 
ram pressure st r ippiii g disc ussed in Secti o n 14.31 This result is consistent with previous studies by 
Catinella et al.l (l2013l ) and IFabello et al.l (120121 ). Nevertheless, the detection ratio does decrease 
significantly in the inner-most regions of even the lowes t mass groups which is cq i isistent with the 
detection of X-rays in some halos of these masses (e.g., iHelsdon fc PonmanI 120001 : iMulchaevIl200fll : 
Zabludoff fc Mulchaev 1998). 


113.35 


Assuming the IGM temperature to be 1 keV for the lowest mass groups {M^ < 10 ^—n ^iviq j 
in this study, our calculation indicates that ram pressure should only be effective for galaxies 
with stellar masses below 10 ® Mq yet all of the galaxies in this sample are more massive than 
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this limit and still show a decreasing detection ratio towards the group cent er. Either the IGM 
temperature in these low m ass groups is lower than assumed (< 1 keV, e.g., iHelsdon fc: Ponman 
2000 : iTrinchieri et al.ll2012l l implying a larger mean IGM density tha n 8 x 10~^cm~^ fro m Equation 


[8] which is similar to the value for poor groups (8.9 x 10 ^cm Zabludoff &: MulchaevI Il998l l 
enough for ram pressure stripping of dwarf galaxies (e.g., iBureau &: CarignanIl2002l l or a process 
other than ram pressu re stripping is responsible for removing the gas in these low mass halos (e.g., 


Rasmussen et al.ll2ni2l ). 


Using a l arge s ample {N ~ 23000) of galaxies with H I mass inferred from optical photometry, 


Zhang et al.l (|2013l l investigate the effect of galaxy stellar mass on the H I content of galaxies in 


groups and clusters and find that lower stellar mass galaxies show a stronger radial decrease in 
H I gas content toward the brightest cluster galaxy and that this trend becomes more significant 
for galaxies with lower stellar surface density. While we can not probe the effect of stellar surface 
density for given galaxy stellar mass due to small sample size, Figure [9] shows that for the lowest 
mass halos it is only the low mass galaxies that are significantly stripped while in the high mass 
halos all of the galaxies are likely to be i mpacted. Our result that the halo mass impacts the degree 
of gas depletion contrasts with those of IZhang et al.l (j2013l ) who find no difference as a function 
of group velocity dispersion as a proxy for halo mass. This dif ference may be related to the large 
dispersion in the velocity dispersion - halo mass relation le.g.. IWeinmann et al.l 1201011 and/ or the 
additional impact of tidal stripping on the gas in their sample galaxies (jZhang et al.ll2013l ). For 
this study galaxies with neighbors within 1.75 ' (which means no pair galaxies within 20 kpc at 
z = 0.01 and 105 kpc at z = 0.055) have been excluded. Many of the excluded gala xies with small 


veloci ty separations are likely to be undergoing tidal interactions as suggested by IScudder et al 


(|2012). Therefore , this s ample should be significantly less strongly impacted by tidal interactions 


than 


these galaxies. 


Zhang et al.l ()2013l ) and should more clearly show the impact of ram pressure stripping on 


5.2. Effect of Pre-processing 


In hierarchical models, the gas content in galaxies d ecreases during interactions within small 
groups before they enter a large halo (e.g., iFuiital 12004 1. Observing this pre-processing of small 
groups is challenging because the current merger rate is considerably lower than in the past 


(jGottlober et al.M2001h . However, there is supporting evidence for pre-processing revealed, for 


the Coma Supercluster ( 

Clvbulski et al. 

20141 

structures of galaxy clusters ( 

MahaiaiJ 

20131: 

of galaxies ( 

Cortese et al. 

20061 '). 


and the ionized gas distribution in compact groups 


The detection ratio of galaxies at the largest group-centric distances (i.e., beyond the virial 
radius) may provide a measurement of how environment has affected the systems prior to their 
falling into the group. The upper left panel of Figure [9] shows with marginal (1.5 cj) statistical 
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significance that the H I detection fraction in massive halos {M^ > /i“^Mq) is 68% of the 

detection fraction for the control sample even beyond the virial radius (i.e., d/Ri^Q ~ 1). The effect 
is even more strongly pronounced for the lowest mass galaxies as shown in the lower right panel of 
Figure [9l where the detection fraction of the lowest mass galaxies is less than that of the control 
sample by 43% corresponding to a 2a difference. 


While not definitive because of the low statistical significance, this depression in the detection 
fraction for small galaxies in massive halos provides a tantalizing indication that these sm all galaxies 


are ak eady H I deficient when they fall into the massive halos, as also suggested by iJaffe et al 
(|2012l i based on the H I measurement of compact group. This trend of gas deficiency at the 
outskirt of massive hal os supporting pre-p rocessing scenario is predicted by semi-analytic model 


Rafieferantsoa et al 


of mass accretion fe.g..lMcGee et al.N2nn9l ) and hydrodynamic simulation of H I gas content (e.g., 


201411 ■ and also iinplied by recent observational studies using statistical sample 


fe.g.. lHess fc WilcotsI 2013 : Hou et al. 2014) 


In considering the significance of this depression at large group centric distance and the other 
results presented here, it is important to consider whether the sample selection criteria (excluding 
groups with fewer than 5 members, excluding galaxies with close neighbors, and excluding galaxies 
with fewer than 8 controls) might bias these results. While bias can not be ruled out, these selection 
criteria should not influence these results within the bounds of the sample (i.e., we can not evaluate 
the impact of the group environment on the smallest and most compact groups) with the exceptions 
of the close neighbors. The removal of galaxies with close neighbors will preferentially eliminate 
interacting pairs. Interacting pairs, however, have not shown any significant decrease or increase 
in their detection fraction (Fertig et al. 2015 submitted) implying that the impact of removing 
galaxies with close neighbors will not introduce a significant bias to the detection ratio estimate. 


6. CONCLUSIONS 

We have investigated the H I gas content in groups and clusters over a wide range of halo 
masses and the results of this study are: 

• There is no statistically significant variation in H I gas-to-stellar mass ratio with projected 
distance from the group center in the sample studied here. 

• The detection ratio of galaxies decreases towards the center of groups spanning the range of 
masses studied here with the most significant decrease occurring towards the centers of the 
most massive halos. 

• We frame the GG criteria in terms of halo mass and stellar mass and find that the change in 
the radial distribution of detection ratio with respect to halo mass is qualitatively consistent 
with ram pressure stripping while the decrease in the detection ratio for galaxies with M* < 
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10 ®'® Mq in the lowest mass halos is not consistent with ram pressure stripping unless the 
temperature of the IGM is lower than assumed (i.e., < 1 keV). 

• The detection ratio for low (M* < 10®'®Mq) stellar mass galaxies just beyond the virial radius 
of the most massive groups is 2a below unity, which could be due to pre-processing of these 
galaxies before they fall into the larger groups. 


The lack of a significant detection of H 1 depletion in groups and clusters, indicates that AL¬ 
FALFA is not sufficiently deep to measure this effect because of the large galaxy to galaxy dispersion 
in gas content. Even deeper H 1 surveys like the Arecibo Galaxy Environment Surve y (AGES) do 


not f ind strong evidence for H 1 deficiency inc reasing toward the center of the Virgo (iTavlor et al 


2013) and Goma/Al367 (jGortese et al.l 120081 ) clusters. Explicitly measuring the amount of envi¬ 


ronment dependent gas depletion wi ll require significantly deeper large statistical samples which 
might have to wait for the full SKA ( Gatinella et al.ll2013l ). 
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A. A simple form of the Gunn and Gott criterion depending on halo mass 


A simple criterion for ram pressure stripping first proposed by I Gunn fc Gotti (jl972l l is 


PicmV > 27rGS*I]isM 


(Al) 


where /?icm is the density of hot intracluster medium (IGM), v is the speed of infalling galaxy, and 
S* and SjsM is the density of galaxy stellar mass and interstellar medium (ISM) respectively. Using 
a simple assumption that the galaxy stellar mass (M*) and ISM mass (Mism) are distributed along 
the same characteristic radius rg, we can rewrite this criterion as 

( _ 1 _ ] ( _^_ 

VlO ^cm 3/ \ 1000km s ^ 

where n is the mean electron number density in a galaxy cluster such that: picM = nnip and a is 
the line of sight velocity dispersion of galaxies in groups and clusters (i.e., 3 cj^). 


2 / 
> 0.37 ^ 


iylOiOMc 


M, 


ISM 


M* 


lOkpc 


-4 


(A2) 


If one defines the halo mass = Migo and the halo virial radius Rmr = Riso using the 
quantities at an overdensity of 180 and combines the virial theorem cr^ = and Equation [3l 

_ -txvir 

one can modify Equation IA2I and the ram pressure criterion including halo mass is written as 


/ ^ \ / Mh 

\10“3cm“3/ yiO^^h'^M© 


2/3 

>4.23 


M* 

lOiOMo 




(A3) 


where we replace Mjsm by cold atomic gas mass Mhi. 


Eurthermore, if we adopt the scaling relations between M* and Vg and between M* and Mhi, 

the physics on the right hand side of the above criterion can be described by galaxy stellar mass 

/ \ - 0.86 

AT* only. Among galaxy scaling relations, we use the relation Mhi/M* Ri 0.72 ( ) 


Eabello et al 


2011 


from [Masters et al.l (l2012li which finds the similar relation to other studies (|Catinella et al. 


201C 


Toribio et al.l 120111 1 and the relation between M* and Vg derived from the re¬ 


lation between galaxy stellar mass and half-light radius for late type galaxies with Sersic index 

0.14 / \ 0.25 


- < 2-5, 


0-1 te) 


1 + 


M, 


3.98xl0i“M© 


m 


Shen et al.l ( 20031 1. We define as 5 x rso 


which encloses nearly 100 % and 90% of t he total luminosity from the Sersic profile for re = 1 and 
re = 4 respectively (e.g., lYoon et al.l 120111 ). Then we can simplify the Equation IA3I as follows. 


re 


10 3cm 3 


Mh 


IO14/1-IM0 


2/3 


> 1.24 


M* 


(lOiOMo 


0.58 


1-F0.25 X 


M* 


IOIOMq 


-1 


(A4) 


On the other hand, soft X-ray emission from galaxy cluster can be used to infer the mean 
electron number density re of the cluster. One of the robust scaling relations is the Yx-M^qq 
relation (jArnaud et al.l 120071 ): 





Mg^50oTx 

W^^M^keV 


0.548 


(A5) 
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where h?‘{z) = Tx is the kinetic temperature of ICM, M 500 is the total mass 

of a cluster enclosed by the radius -R 500 and Mg^soo is the mass of hot gas enclosed by -Rsoo- By 
substituting Mg^soo with ^R^Q^ump, we express the n given halo mass and temperature at z = 0 
(i.e., h{z) = 1) as follows 


n 


10 ^cm 3 


o.esh^ 


T f Miso 

keVy 


0 


0.824 


(A6) 


wher e Migo is converted from M 500 (IReiorich et al.ll2013li by assuming NFW profile (INavarro et al 
19971 1 . By combining Equation IA4I and IA6I and using ~ Migo, one can find a region where ram 


pressure stripping becomes efficient for given galaxy stellar mass, halo mass and kinetic temperature 
of ICM in the halo, as shown in Figure [TTJ 
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Fig. 1.— Distribution of the galaxy sample from cross-matching between the NASA Sloan Atlas, 
the SDSS group and the a.40 catalog in the H I flux and line width plane. Triangles (blue in the 
online version) are the group samples and grey dots are the contro l samp les. Lines (red in the 
online version) are completeness limits determined by iHavnes et al.l (1201 il l. For both the group 
sample and the control sample approximately 41% and 6% of galaxies are under the 90% and 25% 
completeness limits respectively. 
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Fig. 2.— Stacked position of group member galaxies relative to the center, normalized by each 
group halo’s i?i 80 - The grey dots are the SDSS galaxies within the a.40 catalog footprint and 
triangles (blue in the online version) are the ALFALFA detected galaxies. The dotted line represents 
i?i8o- Left: galaxies in groups with Mq. Right: galaxies in groups with > 

IQ^smh-iMQ 
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Fig. 3.— Image of AGC8842 with Hickson Compact Group 69 galaxies within 1.75 ' marked with 
arrows. The ALFALFA beam, shown as the blue circle here, includes multiple galaxies which 
boosts the gas-to-stellar mass ratio of the central galaxy identified as the best candidate for the 
ALFALFA source. We avoid such cases by removing groups that have ALFALFA sources that can 
be associated with multiple galaxies within a 1.75 ' radius. 
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Fig. 4.— Standard deviation of the average log(MHi/M*) for the control samples for each pair of 
the 50 stellar mass and redshift values, as function of the number of control samples. The different 
black lines are 20 separate calculations of the above values for standard deviation as a function of 
control sample size. The grey (red in the online version) line is the ‘average’ of these 20 lines. The 
thin solid (red in the online version) line has a slope of the statistically expected change 

in the standard deviation as the number of controls change. The line has been placed with an 
arbitrary offset to show that the data matches the expected slope well for samples with 7 or more 
controls. However, for large numbers of controls the number of groups that have a full complement 
of controls drops rapidly. We find that the optimal control sample size is 8. 
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Fig. 5.— Stellar mass and redshift for galaxies in groups represented by circles (blue in the online 
version) and for isolated galaxies represented by small black points. The filled circles (red in the 
online version) are the galaxies in the final sample that have a full complement of 8 controls. 
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Fig. 6.— The difference between logM* for the galaxies in the group sample and the log of 
the average M* for their 8 associated control galaxies as a function of the projected group centric 
distance for different halo mass bins. Projected distance is normalized by group virial radius (i?i8o)- 
Data points (magenta in the online version) are individual galaxies and the associated error bar is 
the range of the value determined by the 8 control samples. The solid line represents a mean of 
the AlogM* in each radial bin. The dashed lines enclose the 5*^ and 95*^ quantile of the radial 
distribution generated by random shuffling of galaxy radial position. Note that no significant radial 
variation of the stellar mass relative to the controls (AlogM* < 0.01) is seen for galaxies in all 
halo mass bins. 
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Fig. 7.— Same as Figure [6] for the difference of redshift between the galaxies in the group sample 
and the mean of the corresponding 8 control galaxies, as a function of the projected group centric 
distance for different halo mass bins. As in Figure [6l no significant radial variation of the redshift 
relative to the controls {Az < 0.0002) is seen for galaxies in all halo mass bins. 
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Fig. 8.— The difference in H I gas-to-stellar mass ratio between galaxies in the group sample 
and the mean of the corresponding 8 control galaxies, as a function of the projected group centric 
distance for different halo mass bins. Symbols and lines are the same as those in Figure [H 
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Fig. 9.— Upper left panel: The variation of the detection ratio as a function of the projected group 


centric distance for different halo mass bins. Other panels: The detection ratio for each halo mass 


bin with the detection ratio for two groups of galaxies divided by galaxy stellar mass. Open points 
connected by red and blue lines show these relationships for massive galaxies (M* > IO^'^Mq) and 
low mass galaxies (M* < Error bars are Poisson errors. The detection ratio decreases 

toward group center and the decrement is most dramatic in large mass halos. 
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Fig. 10.— Upper left panel: The variation of the detection ratio as a function of the projected 
group centric distance for different halo mass bins. Other panels: The detection ratio for each halo 
mass bin for galaxies separated by galaxy Sersic index. Open points connected by red and blue lines 
show these relationships for early type galaxies (n > 2) and late type galaxies (n < 2). Control 
samples for each early and late type galaxy are selected from the pool of isolated field galaxies with 
n > 2 and n <2 respectively. Error bars are Poisson errors. The detection ratio decreases toward 
group center for both early and late type galaxies. 
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Fig. 11.— The halo mass-galaxy stellar mass plane that determines the conditions for efficient 
ram pressure stripping. Black lines represent different ICM temperatures based on Equation [7] 
and divide the halo mass-galaxy stellar mass parameter space into the regions where ram pressure 
dominates and where galaxy gravity dominates. Ram pressure stripping is efficient in the region 
above the black lines for a given ICM temperature. The distribution of the group galaxy sample 
used in this study is shown by black points, the colored squares represent the weighted average of 
Alog(MHi/M*) in each bin. 
















